INVERSE SCATTERING TRANSFORM FOR THE TODA 
HIERARCHY WITH QUASI-PERIODIC BACKGROUND 
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Abstract. We provide a rigorous treatment of the inverse scattering trans- 
form for the entire Toda hierarchy in the case of a quasi-periodic finite-gap 
background solution. 



1. Introduction 

Since the seminal work of Gardner et al. [7] in 1967 the inverse scattering trans- 
form is one of the main tools for solving completely integrable wave equations and 
numerous articles have been devoted to this subject since then. However, with 
only a few exceptions, only the case of short range perturbations of a constant 
background solution has been investigated. The first to consider a non-constant 
background seem to be Kuznetsov and A.V. Mikhailov, [H], who informally treat 
the Korteweg-de Vries equation with the Weierstrafi elliptic function as stationary 
background solution. Formulas for arbitrary periodic background were given by 
Firsova in [5]- However, the discontinuous behavior at the band edges is not inves- 
tigated there and neither is the question when the time dependent scattering data 
satisfy the hypothesis necessary for the Gel'fand-Levitan-Marchenko theory (see, 
e.g., JOI f° r the constant background case). Moreover, the normalization chosen 
there blows up when the eigenvalues lie at certain positions in the gaps. 

Our aim in the present paper is to provide a rigorous treatment of the inverse 
scattering transform for the Toda hierarchy in the case of quasi-periodic finite- 
gap background solutions (including the case of arbitrary periodic solutions). We 
include a discussion of the problems arising from the poles of the Baker- Akhiczcr 
functions at the Dirichlet eigenvalues (which are absent in the constant background 
case). In the constant background case this was first done for the Toda equation 
by Flaschka in |Sj on an informal level. The first rigorous treatment for the entire 
Toda hierarchy was given by one of us in ^2] ( see a ls° EH- |15| . |14j1. For further 
results on solutions of the Toda lattice with nontrivial spatial asymptotics see pQ. 

After introducing the Toda hierarchy in Section [21 we will first show that a 
solution will stay close to a given background solution in Section [3] This result 
implies that a short range perturbation of a quasi-periodic finite-gap solution will 
stay short range for all time and it shows that the time dependent scattering data 
satisfy the hypothesis necessary for the Gel'fand-Levitan-Marchenko theory. This 
result constitutes the main technical ingredient for the inverse scattering transform. 
In Section 01 we review some necessary facts on quasi-periodic finite-gap solutions 
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and in Section0we compute the time dependence of the scattering data and discuss 
its dynamics. Finally we establish the conserved quantities of the Toda hierarchy 
in this context. 

2. The Toda hierarchy 

In this section we introduce the Toda hierarchy using the standard Lax formalism 
(OH)- We first review some basic facts from [2] (see also [T4*]L 
We will only consider bounded solutions and hence require 

Hypothesis H. 2.1. Suppose a(t), b(t) satisfy 

a(t) G £°°(Z, E), b(t) G £°°(Z, E), a(n, t) ^ 0, (n, (JeZxR, 

and let 1 1-> (a(t),b(t)) be differentiable in £°°(Z) £°°(Z). 

Associated with a(t), b(t) is a Jacobi operator 

(2.1) H(t):£ 2 (Z)^£ 2 (Z), f~T(t)f, 
where 

(2.2) T(t)f(n) = a(n, t)f(n + 1) + a(n - 1, t)f(n - 1) + 6(n, *)/(n) 

and i? 2 (Z) denotes the Hilbert space of square summable (complex- valued) sequences 
over Z. Moreover, choose constants Co = 1, Cj, 1 < j < r, c r+ i = 0, set 

j 

3j(M) = ^Cj^i{5 n ,H{t) l 5 n ), 

1=0 

i 

(2.3) hj(n,t) = 2a(n, t) Cj-e(5 n+1 , H(t) e 5 n ) + c j+1 , 

£=o 

and consider the Lax operator 

r 

(2.4) P 2r+2 (t) = -H{ty +1 +J2(Mt)9j(t)S+ - hjityHity-S+gr+itf), 

3=0 

where S' ± /(n) = f(n± 1). Restricting to the two-dimensional nullspace Ker(r(i) — 
z), z G C, of r(t) — z, we have the following representation of P 2r+ 2{t) 

(2.5) P2r+2(t) =2a{t)G r (z,t)S + - H r+1 (z,t), 

Kcr(-r(t)-z) 

where G r (z, n, t) and iJ r +i(z, n, i) are monic polynomials in z of the type 

r 

G r (z,n,t) = ^2z J g r -j(n,t), 

3=0 

r 

(2.6) H r+ i(z,n,t) = z r+1 + ^ z J h r -j(n,t) - g r+ i(n,t). 

3=0 

A straightforward computation shows that the Lax equation 

(2.7) j t H{t) - [JV+aW, #(*)] =0, < G K, 
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is equivalent to 

TL r (a(i), b{t))\ = d{t) - a{t)(g+ +1 (t) - g r+ i(t)^ = 0, 

(2.8) TL r (a(i), b{t)) 2 = b(t) - (h r+1 (t) - = °' 

where the dot denotes a derivative with respect to t and f^(n) = f(n± 1). Varying 
r e N yields the Toda hierarchy Th r (a, b) = (TL r (a, 6)i,TL r (a, b) 2 ) = 0. 

Finally, we recall that the Lax equation ((2.7(1 implies existence of a unitary 
propagator U r (t,s) such that the family of operators H(t), f £ 1, are unitarily 
equivalent, H{t) = U r (t, s)H(s)U r (s,t). 

3. The initial value problem 

First of all we recall the basic existence and uniqueness theorem for the Toda 
hierarchy (see, e.g., [H], [13], or d Section 12.2]). 

Theorem 3.1. Suppose (a ,b ) G M = £°°(Z) © £°°(Z). Then there exists a 
unique integral curve t t— > (a(t),b(t)) in C°°(M, M) of the Toda equations, that is, 
TL r (a(t),b(t)) = 0, such that (a(0),6(0)) = (a ,b ). 

In |12| it was shown that solutions which are asymptotically close to the constant 
solution at the initial time stay close for all time. Our first aim is to extend this 
result to include perturbations of quasi-periodic finite-gap solutions. In fact, we 
will even be a bit more general. 

Lemma 3.2. Suppose a(n,t), b(n,t) and a(n,t), b(n,t) are two arbitrary bounded 
solutions of the Toda system satisfying \S. 1\) for one to <E K, then 18. 1)) holds for 
all f £ R, that is, 

(3.1) ^2w(n)(\a(n,t) - a(n,t)\ + \b(n,t) - b(n,t)\j < oo, 
where w(n) > 0. 

Proof. Without loss of generality we assume that to = 0. Consider the expression 
\\(a(t),b(t))\\* = 5^tfl(n)(|o(n,t) - o(n,t)| + \b(n,t) - b(n,t)\) 

which remains finite at least for small t since there is a local solution of the Toda 
system with respect to ||.||*. We claim that the following estimate holds 

Y,w(n)\g r (n,t)-g r (n,t)\ <C r \\(a(t),b(t))\U, 

(3.2) ^2w(n)\h r (n,t)-h r (n,t)\ < C r \\(a(t), 6(t))||*, 

where C r = C(\\H(0)\\,\\H(0)\\) is a positive constant due to HO Let us prove 
(13.21) by induction on r. It suffices to consider the case where Cj = 0, 1 < j < r, 
since all involved sums are finite. In this case (141 Lemma 6.4] shows that gj(n,t), 
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hj(n,t) can be recursively computed from go(n,t) = 1, ho{n,t) = via 

g j+1 (n,t) = ~(hj(n,t) + hj(n - l,t)) + b(n,t)gj(n,t), 

j ^ j 

h j+1 (n,t) = 2a(n,t) 2 ^2g J _i(n,t)gi(n + l,t) - ^y^hj-i(n,t)hi(n, t). 

1=0 1=0 

Hence 

<7i(n, t) — gi(n, t) — b(n, t) — b(n, t), hi(n, t) — hi(n, t) — 2(a(n, t) 2 — a{n, t) 2 ), 
and 

9r+i(n,t) -g r+1 (n,t) = (b(n,t) -b(n,t))g r (n,t) + b(n,t)(g r (n,t) -g r (n,t)) 
+ — (h r (n, t) — h r (n, t) + h r (n — 1, t) — h r (n — 1, i)) , 
1 r 

h r+1 (n,t) - h r+1 (n,t) = - (h r -i(n,t)hi(n,t) - h r -i(n,t)hi(n,t)) 

1=0 

r r 

+ 2a(n,t) 2 ^g r ^i(n,t)gi(n + l,t)- 2a(n,t) 2 ^g r -i(n,t)gi(n + 

(=0 1=0 

Since both sequences a, b and a, b are solutions of the Toda system, (|2.8|) yields 
a(t) ~ a(t) = a(t)(g+ +1 (t) ~ g r+1 (t)) - a(t)(g+ +1 (t) - g r+1 {t)) 

= Ht) - a(t))(g+ +1 (t) - g r+ i(t)) ~ a(t)(g+ +1 (t) - g+ +1 (t) + g r+1 (t) - g r+1 (t)) 
and thus 

\a(n, t) — o(n, t)\ < \a(n, 0) — a(n, 0) I + C\ I \a{n 1 s) — a(n, s)\ds 



+ \\H(0)\\ / (\g r+1 (n + l,s) - g r+1 (n + l,s)\ + \g r+ i(n,s) - g r+1 (n,s)\)ds, 
Jo 

\b(n, t) - b(n, t)\ < \b(n, 0) - b(n, 0)| 

+ / (\h r +i(n,s)-h r+1 (n,s)\ + \h r+ i(n-l,s)-h r +i{n-l,s)\)ds. 



By (E2J, ||(o(t),6(t))||, < ||(o(0),6(0))||. + C /„* || (a(s), &(*))||.*», where C = (d + 
2\\H(0)\\ +2)C r +i. Finally, applying Gronwall's inequality 

\\(a(t),b(t))\U < ||(a(0),6(0))||.exp (Ct) 

finishes the proof. □ 

4. Quasi-periodic finite-gap solutions 

As a preparation for our next section we first need to recall some facts on quasi- 
periodic finite-gap solutions (again see [5] or [T^)). 

Let M be the Riemann surface associated with the following function 

29+1 

(4.1) R\' g 2 +2 ^) = - II \f z ~ E h E <E 1 <---< E 2g+1 , 

j=o 

where g G N and yC is the standard root with branch cut along (— oo,0). M is a 
compact, hyperelliptic Riemann surface of genus g. A point on M is denoted by 
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p = (z, ±i? 2 ^ 2 ( z )) = (z, ±), z£ C, or p = oo±, and the projection onto C U {00} 
by 7r(p) = z. 

Now pick 5 numbers (the Dirichlet eigenvalues) 

(4-2) {N)U = ^^)U 

whose projections lie in the spectral gaps, that is, fij £ [E2j—i,E2j]. Associated 
with these numbers is the divisor Dp which is one at the points jlj and zero else. 
Using this divisor we introduce 

z{p, n, t) = A po (p) - a po (Dp) - ni oc _ (oo + ) + tU s - E po e C 9 , 
(4.3) z(n, t) = z(oo+, n, t), 

where S po is the vector of Riemann constants, U_ s the ^-periods of the Abelian 
differential fl s defined below, and A po (a po ) is Abel's map (for divisors). The 

hat indicates that we regard it as a (single-valued) map from M (the fundamental 
polygon associated with M) to C 9 . We recall that the function 0(z{p, n, t)) has 
precisely g zeros jlj(n,t) (with /tj(0,0) = fij), where 9{z_) is the Riemann theta 
function of M. 

Then our solution is given by 

2 . 2 0{z{n + l,t))6(z(n-l,t)) 



a q (n, t) = a 



6(z(n,t)) 2 

, . .s , . ~ 1 . d _ / 9(w + z(n.t)) 

(4.4) ^(n,*) = & + 5>(s)— ln( - - 



_ dwj \6(w + z(n — l,t)) 



w=0 



The constants a, b, Cj(g) depend only on the Riemann surface (see ^] Section 9.2]). 
Introduce 

<Pq{P,n,t) = G(n,t) — — -— exp I / uj^ 00 _), 

(4.5) ^,(p, n, t) = C(n, 0, t) ^|^|] exp (n [ Woo+i00 _ + t [ 0.) 



* h 

Po J P0 



where C(n,t), C(n,0,t) are real-valued, 

(4.6) C(M) 2 = ! fe S n ' M !!, C(n,0,t) 2 - W°))^(-1,0)) 



(z(n + l,t))' v > . 7 e{z{n,t))9{z{n-l,t)Y 

and the sign of C(n,t) is opposite to that of a q (n,t). u) 00+j00 _ is the Abelian 
differential of the third kind with poles at oo + respectively 00 _ and Q s is an Abelian 
differential of the second kind with poles at oo + respectively oo_ whose Laurent 
expansion is given by the coefficients (j + l)c s -j associated with TL S (see |14l 
Sections 13.1, 13.2]). Then 



Tq(t)ip q (p,n,t) = n(p)ip q {p,n,t) 
—ip q (p,n,t) = 2a q (n,t)G s (p, 
(4.7) = P q , 2s+ 2(t)i>q(p,n,t)(n) 



—ipq(p,n,t) = 2a q (n,t)G s (p,n,t)ip q (p,n + l,t) - H s+1 (p,n,t)ip q (p,n,t) 
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It is well known that the spectrum of H q (t) is time independent and consists of 
g + 1 bands 

9 

(4-8) a{H q )={J[E 23 ,E 2]+1 ]. 

For further information and proofs we refer to ^] Chapter 9] . 
Next let us renormalize the Baker- Akhiezer function 

such that ipq(p, 0, t) = 1 and let us abbreviate 

(4.9) exp (a s (p, t)) = 0, t) = C(0, 0, t) ^'Xll ex P (* 

One can show ([H], (13.49)) 



(4.10) 
where 



exp (a s (p,t)) =exp|^y (2a 9 (0, x)G s (p, 0,^)^(^,0, x) - ff s+ i(p,0, a;)) dx 
'G fl (p,0,t) / 1/2 f 4 G.Cp.O.x) ' 

ex P ^ 2g+2 (P)y o Gg(P;0;X) ^ 



G ff (p,0,0) 



(4.11) G ff (2,n,t) = IJ(2!-/i J -(n,t)). 

3=1 

However, let us emphasize that the integrals in (|4.1(J|I are only well-defined until we 
hit the first pole of the integrand. After that time they have to be interpreted as 
principal values (compare also Remark [O] (i)), respectively (|4.9|) has to be used. 
In addition, observe 

ot s ,+ {z,t) + a Sj _(z,t)J = - ^ - — , 



(4.12) ex P (« s , + (M)-« 8 ,-fet)) =expf2ii 2 %(z) T 



dx I . 



For A 6 cr(Hq) we have 



(4.13) a.,±(A,t) =a» )T (A,t), 

but note that a SiT (A,t) 7^ T a s.±(A, £) unless 5 = 0. 

5. Inverse scattering transform 
Let a(n,t), b(n,t) be a solution of the Toda system satisfying 

(5.1) M (H n ' *) - a «( n ' *)l + *) - b q (n, i)|) < 00 

for one (and hence for any) to G M. In (see also JB]) we develop scattering theory 
for the Jacobi operator H associated with a(n), b(n). Jost solutions, transmission 
and reflection coefficients depend now on an additional parameter i £ R, More 
precisely, 

(5.2) a{H{t)) = a{H), a ess (H) = a(H q ), a p (H) = { Pj }U C M\a(H q ), 
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where q G N is finite. The Jost solutions ip±(z, n, t) are normalized such that 

ip±(z, n, t) = ip q> ±(z, n, t) (1 + o(l)) as n — > ±00. 

Transmission T(X, t) and reflection R±(\, t) coefficients are defined via the scat- 
tering relations 

(5.3) T(\,t)4>±(\,n,t) = 4> T (\,n,t) + R T (\,t)ip T (\,n,t), A G cr(H q ), 
and the norming constants 7± t j(t) corresponding to pj G a p (H) are given by 

(5.4) 7± J (r 1 = Ei^> n > t )i 2 - 

To avoid the poles of the Baker-Akhiezer function, we will assume that none of 
the eigenvalues pj coincides with a Dirichlet eigenvalue jtXfc(0, 0). This can be done 
without loss of generality by shifting the initial time to = if necessary. 

Remark 5.1. Due to this assumption there is no need to remove these poles for the 
definition of r y±j, as we did in Since the Dirichlet eigenvalues rotate in their 
gap, the factor needed to remove the poles would only unnecessarily complicate the 
time evolution of the norming constants. Moreover, these factors would eventually 
cancel in the Gel'fand-Levitan-Marchenko equation, which is the only interesting 
object from the inverse spectral point of view in the first place. See Remark \5.4\ 
below. 

Lemma 5.2. Let (a(t),b(t)) be a solution of the Toda hierarchy such that \5.1)) 
holds. The functions 

(5.5) ip±(z,n,t) = exp(a St ±(z,t))ip±(z,n,t) 
satisfy 

(5.6) H(t)ijj±(z,n,t) = zijj±(z,n,t), ^-tp ± (z,n,t) = P 2s+ 2(t)ip±(z,n,t). 

dt 

Proof. We proceed as in | 12l Theorem 3.2]. The Jost solutions ip±(z,n,t) are con- 
tinuously differentiable with respect to t by the same arguments as for z (compare 
; 3,i Theorem 4.2]) and the derivatives are equal to the derivatives of the Baker- 
Akhiezer functions as n — > ±00. 

For z G p(H), the solution u±(z, n, t) of (|5.()|) with initial condition ip±{z 1 n, 0) G 
£j_(Z) remains square summable near ±00 for all t G K (see ^21 or [HI Lemma 
12.16]), that is, u±(z, n, t) = C±(t)ip±(z, n, t). Letting n — * ±00 we see C±(t) = 1. 
The general result for all z G C now follows from continuity. □ 

This implies 

Theorem 5.3. Let (a(t),b(t)) be a solution of the Toda hierarchy such that \5.1)) 
holds. The time evolution for the scattering data is given by 

T(\,t) = T(A,0), 
R±(X,t) = fl±(A,0)exp(±(a s , + OM)-a s ,-OM)))> 
(5-7) 7 ±J -(t) = 7 ±ii (0)exp(2a s , ± ( / ; i ,t)), 1 < j < q. 
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Proof. Since a Sj ±(A,i) = a SiT (X,t) by (|4.13|) . we have 



W(^(X,t)^+(X,t)) 



exp(g a ,_ (A, t)) exp(a,, + (A, t)) W{$- (A, t), g_ (A, t)) = 
exp(a s ,_(A,t))expK_(A ) i)) W(^-(A,i),V+(A,t)) [ 1 >' 

Here we used that the Wronskian W n (f,g) = a(n){f(n)g(n + 1) — /(n + l)g(n)) 
of two solutions satisfying l|5.6|) does not depend on n or t (see Eli EI Lemma 
12.15]). Similarly, 



exp(q a ,±(A, t)) exp(a SiT (A, t)) W/Q T (A, t), ip± (A, t)) 
± ' exp(a s , T (A, t)) exp(a s ,±(A, f )) W{$± (A, t), V T (A, t)) 

= exp(a a ,±(jz,t) - a StT (z,t))R±(X,0). 

The time dependence of j± j(t) follows from || ip±(pj, ., £) | = ||£/ s (i, 0)^>±(/9j, ., 0)|| = 
||^±(p i ,.,0)||. ' ' □ 

A few remarks concerning this theorem are in order: 

Remark 5.4. (i) First of all note that the time dependence of R±(X,t) is not 
continuous at the band edges (except the lowest and highest). In fact, the integrand 
in \4-.l'^jj has a pole whenever pj(0,t) hits a band edge. Suppose that pj(0,t) gets 
close to £?2j att = t\, then (assuming G a (E2j,0,ti) ^ for simplicity) 



nj(p,t) = e 2j + vj (t - t x y + o(t - h) 6 , 



Ili^O*j(o,ti)-w(o,*i)) 



As long as t <t\ we have 



lim (a» )+ (A,t)-a.,_(A,i)) = 0. 



But at t = t\ we obtain (setting X = E 2 j + s 2 and using v.j = (— I) 9 ^\\vj\) 
A ligK + (A, «)-»„- (A, t)) - to M&Wjf 



lim 2 / 1 -= ^ ^dx = 2(-l) 9 - J 'ilimarctan(^i) = (-1) 



e 



T/iis is m agreement with the conditions 

lim R 1/2 r\ ^±(V)+ 
Um K 2r+2 (AJ — 

lim i?2, 7 + 2 (A) ^g^}" 1 = °> & = H (°> *) 



lim i^ 2 (A) =0, £?, # Mj (0, t) 



( 5 ' 8 ) ,• nV2 / A x J?.±(A,t)-l _ 
™ ^2r+2( A ) T(I) 



found in 

fn) XTie norming constants j± j(t) will vanish whenever (pj,±) = pk{0,t). TTizs 
is important, since the norming constants must cancel the pole of ip q ,±(pj,n,t) 
at these points such that the Gel'fand-Levitan-Marchenko equation remains well- 
defined. 

In fact, these poles come from the normalization ip q ,±{z,0,t) = 1 and this nor- 
malization is natural in the constant background case since it renders ip q _±(z,n,t) 
time-independent. In the ( quasi- )periodic background case this is no longer true 
and it seems more natural to allow a time- dependent normalization and replace 
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il>q,±{z, n -, t) by the time- dependent Baker- Akhiezer functions ip q ,±(z, n, t) rendering 
the norming constants and the reflection coefficient time-independent. 

(Hi) Finally let us discuss the connection with the different normalizations for 
j± t j used in |3j respectively 

In [3] we have multiplied ip q ,±(X, n, t) by a factor to cancel the poles at the Dirich- 
let eigenvalues. However, since the Dirichlet eigenvalues change sheets during time 
evolution, this normalization renders tp qi ±(X,n,t) discontinuous with respect to t 
and complicates the time evolution of the corresponding norming constants: 

(5.9) j ±d (t) = 7 ±li (0) exp ( ± 2RlH 2 ( Pj ) J* ^El^dx 



± y t* 2GMx),Q,x)R^; +2 (Mx)) ^ 
fzw o (Pi ~ M 0*0) flft^ (« ( x ) ~ Mfc ( x )) 



This clearly justifies our decision to use the freedom of choosing the initial time to 
avoid the poles in the first place. 

If H q is periodic with period N , the normalization chosen in |S] is to divide 
ip q .±(X,n,t) by the square root of jj^2n =1 ip q ,+ (z,n)ip q ^_(z,n). To get a well- 
defined real-valued result this object should be positive. However, since one can 
compute (see ^Jj\) 

N N-l . 

(5.10) ^>+( 2 ' n )%,-(*, n) = N Y[ Z —±, 

ra=l j=l Z M 

where Xj € [E^j—i, E%j] are the zeroes of uj 00 __ 00+ , this factor can be negative or 
even zero. While this choice does give a simple time dependence 

(5.11) J±,j(t) = l±,j{0) exp(a Si± (pj,t) - a s ^(pj,t)), 

it does not cancel the poles at the Dirichlet eigenvalues ( it just renders poles into 
square root singularities) plus it introduces additional singularities at Xj. In par- 
ticular, it is unsuitable if one of the eigenvalues coincides with some Xj . 

In summary, since Lemma 13 . 21 ensures that (I5.1|) remains valid for all t once it holds 
for the initial condition, we can compute i?±(A, 0) and 7±.j(0) from (a(n, 0), 6(n, 0)) 
and then solve the Gel'fand-Levitan-Marchenko (GLM) equation to obtain (a(n, t), b(n, t)) 
as in In this respect the most important ingredient is the time dependence of 
the kernel of the GLM equation which follows from Theorem 15.31 

Theorem 5.5. The time dependence of the kernel of the Gel'fand-Levitan-Marchenko 
equation is given by 

F ± {m,n,t) = -Rc / R± (A, 0)ip q ,± (A, m, t)ip q ± (A, n, t) Gg ^ °' 0) dX 

q 

(5.12) + Y 7± ,3 (0)i>q,±{pj,m, t)^q,± {Pj,n, t). 

i=i 

Finally note that since the transmission coefficient is time independent we obtain 
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Corollary 5.6. Let (a(t),b(t)) be a solution of the Toda hierarchy such that \5. 1\) 
holds. Then the quantities 

w A -,kM 

and Tj = tr(/P(t) - H q (t) j ), that is, 

n = ^ b(n,t) - b q (n,t) 

n£Z 

t 2 = 2(a(n,t) 2 - a q {n,t) 2 ) + {b(n,t) 2 - b q (n,t) 2 ) 
etc., are conserved quantities. 

Proof. This follows from the fact that the transmission coefficient is the perturba- 
tion determinant (in the sense of Krein) of the pair H q and H (see ^^)- ^ 
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